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Abstract 

We construct static, asymptotically flat solutions of 5*17(2) Einstein- Yang-Mills theory 
i-C ■ in 4 + 1 dimensions, subject to bi-azimuthal symmetry. Both particle-like and black hole 

solutions are considered for two different sets of boundary conditions in the Yang-Mills sector, 
^s^j , corresponding to multisolitons and soliton-antisoliton pairs. For gravitating multi-soliton 

■ solutions, we find that their mass per unit charge is lower than the mass of the corresponding 

OO ' un it charge, spherically symmetric soliton. 

o ' 

CO ; 1 Introduction 
iri 

The last years have seen an increasing interest in the solutions of Einstein equations involving more 
than four dimensions. The results in the literature indicate that the physics in higher-dimensional 
general relativity is far richer and complex than in the standard four-dimensional theory. 

Naturally, most of the studies in the literature were carried out for vacuum solutions or to 
configurations with an Abelian matter content. At the same time, a number of results in the 
literature clearly indicate that solutions to the Einstein equations coupled to non Abelian matter 
fields possess a much richer structure than in the U(l) case (see [T] for a survey of the situation 
in four dimensions and the more recent review [2] for d > 4), most notably in that they are not 
restricted to black holes, but can also be regular. 

Physically reasonable stationary vacuum solutions in higher dimensional spacetimes, d > 4, 
fall in two categories, distinguished by their asymptotic behaviours. In the first category, there are 
the static spherically symmetric solutions generalising the d — 4 Schwarzschild black hole, found 
by Tangherlini a long time ago [3], the rotating Myers-Perry solution [4] generalising the four 
dimensional Kerr black hole, and more recently the black ring solutions [5 ,6 . In all these cases, the 
e?-dimensional spacetime approaches asymptotically the Ai d Minkowski background. The second 
category are the black string solutions, and the corresponding black p— branes generalizations [7J. 
The black strings approach asymptotically d — 1 dimensional Minkowski-spacetime times a circle, 
x S 1 , and in the simplest case present translational symmetry along the extra-coordinate 
direction. (Such configurations are important if one supposes the existence of extra dimensions in 
the universe, which are likely to be compact and described by a Kaluza-Klein (KK) theory.) 

As is the case with the usual Schwarzschild black hole, all these vacuum solutions can be 
extended to describe configurations with an Abelian matter content. The inclusion of non Abelian 
matter fields is less systematic and is complicated by the fact that all known such solutions can only 
be evaluated numerically, starting from the earliest found Einstein- Yang-Mills (EYM) solution in 
four spacetime dimensions discovered by Bartnik and McKinnon [5] . 
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Spherically symmetric solutions to EYM systems in d-spacetime dimensions, approaching 
asymptotically the M. d Minkowski background, were constructed systematically in [9]-[TH]. The 
Yang-Mills (YM) sector of the systems studied there consisted of all needed terms belonging to the 
YM hierarchy [THl [T7] , which are higher order in the YM curvature in the manner of the Skyrme 
model. Such terms may arise in the low energy effective action of string theory [TH] [19l [20]. It 
has been established that only in the presence of these higher order in the YM curvature terms, 
does the EYM solution lead to a finite mass. In the absence of such Skyrme-like terms, for ex- 
ample in [211 [22] (in d = 5), the mass of the solution diverges. Both particle like and black hole 
solutions were constructed. The properties of these configurations are rather different from the 
familiar Bartnik-McKinnon solutions [8 in d = 4, and are somewhat more akin to the gravitating 
monopole solutions to EYM-Higgs system [53], which is not surprising since the latter features 
the dimensionful vaccuum expectation value (VEV) of the Higgs field, while the former contain 
additional dimensionful terms entering as the couplings of the higher order YM terms. 

As for solutions to the EYM system in d dimensional spacetime whose vacuum has the struc- 
ture of ftA^ 1 x S 1 like the black string solutions, these are only constructed if one of the spacelike 
dimensions is supposed to be compact, and a Kaluza-Klein descent is performed, essentially elim- 
inating that coordinate. Such solutions are given for d — 5 in [21j-|27j. However, in the present 
work we will not be concerned with this type of solutions. 

Our aim in the present work is to extend the construction of asymptotically flat finite mass 
EYM solutions vacuum, relaxing the constraint of spherical symmetry in the d — 1 dimensional 
spacelike subspace as in [9]-|15j. 

The simplest possibility is to consider the imposition of a symmetry which leads to two a di- 
mensional reduced effective system, rather than the one dimensional one in the previous examples. 
This is the first such attempt in the literature, and the numerical work of solving a two dimensional 
EYM boundary value problem is a task of considerable complexity. To achieve a two dimensional 
subsystem, we have found that the simplest option is to impose bi-azimuthal symmetry on the 
d = 5 static EYM system. This is why we have restricted to d = 5, for otherwise a similar appli- 
cation of azimuthal symmetries in each plane would result in multi-azimuthal Q subsystems, with 
higher dimensional boundary value problems to be solved, technically beyond the scope of this 
work. Indeed, as a warmup for the task at hand, we have carried out the same program in |29j 
recently, with the dilaton replacing gravity. 

While we have restricted to five dimensional EYM solutions for technical reasons, this example 
is of considerable physical relevance since it enters all d = 5 gauged supergravities as the basic 
building block and one can expect the basic features of its solutions to be generic. Also special 
about d = 5 gravitating YM is the particular critical properties of the solutions present in all 
d = Ap + 1 analysed in [T^ , and first discovered in [TO] . Indeed in the d — 5 YM-dilaton ( YMd) 
system, studied in [29] . these critical properties were present, providing yet another confirmation 
that dilaton interactions with YM, mimic [30] those with Einstein gravity. 

The purpose of this paper is to present numerical arguments for the existence of a class of static 
d = 5 solutions to the EYM equations of the model studied in [10], but now, subject to bi-azimuthal 
symmetry. These configurations present a spacetime symmetry group R x U(l) x U(l), where R 
denotes time translation symmetry and the C/(l) factors the rotation symmetry in two orthogonal 
planes. We present both regular and black hole solutions. In the particle like case we find solutions 
with many similar properties to those of the four dimensional SU(2) YM multi-instantons and 
composite instanton-antiinstanton bound states with U(l) x U(l) symmetry, reported in Ref. [31 . 
Dilatonic generalizations of these solutions have been considered in [59], in which higher order 
gauge curvature terms were included in the action to enable the existence of finite mass solutions. 

1 If one applied instead, spherical symmetry in the d — 2 dimensional subspace of the d — 1 spacelike dimensions, 
then the residual subsystems will always be two dimensional irrespective of the value of d. For example in d = 5 
this would be the SO (3) symmetry in the 3 dimensional subspace of the 3 dimensional subspace of the 4 spacelike 
dimensions, exactly as for the axially symmetric instantons |28j . While this may appear to be an attractive 
alternative, we have found that tackling the boundary value problem in that case is a considerably harder task, 
even in d = 4. 
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2 The model 



2.1 The ansatz and field equations 

We consider the five dimensional SU(2) EYM action 



s= I <PxV=g [T^-Cm], CD 



16nG 

where C m is given by the superposition of the p — 1 and p — 2 terms in the YM hierarchy [lOj 

£m = 2 . 2!^ 1 ^ ^ v 2 4!^ 1 ' ^ ^ v P a ' ^ 

with = <9 AI „4„ — + i[Afj,, A v ] the 2-form YM curvature and T^ vp(T = {J-^^pa]} the 

4-form YM curvature consisting of the totally antisymmetrised product of two YM 2-form YM 
curvatures (the bracket [vpa] implies cyclic symmetry). T\ and ti are dimensionfull coupling 
strengths. 

Variation of the action fl} with respect to the metric and gauge potential A^ leads to the 
EYM equations 

JV - l^uR = 8^G (TW + T $ ) , (3) 
1 
2 



nV" + ^r 2 {^ pff , D^"?*} = 0, (4) 



where 



7$ = Tr{^(2p) AlAlA2 ... A2p _ 1 ^(2p)/ lA - A2 - 1 - ^(2p) AlA2 ... A2p ^(2p) AlA - A ^}, (5) 

is the energy-momentum tensor for the p— th YM term in p — 1, 2. 

We consider a d = 5 static metric form with two orthogonal commuting rotational Killing 
vectors 

d,s 2 = -f(r, 8)dt 2 + ^(dr 2 + rW) + sin 2 ^ + ?M r 2 cog2 ^2 (g) 

/(»",») /O*,0) 

where r is the radial coordinate, and 9,<p,ip are Hopf coordinates in 5 3 , with < 8 < ir/2 and 
< ip, ip < 2tt. 

The construction of a YM Ansatz compatible with the symmetries of the above line element has 
been discussed at length in [29] , [31] . The purely magnetic gauge connection has six nonvanishing 
components and reads 

A = i-u 3 a r (r, 9)dr + ^u 3 a e (r 7 9)d9 (7) 



/ 1 1 n \ ( 1 1 n 

+ ( ^u lX Hr, 9) + 2 U 2X 2 (r, 8) + -u 3 dip + f -«i^(r, 0) + -u 2 f 2 (r, 0) + 7^3 



where ui = sm.n((p + i/>)ci — cos n{ip + ip)c2, ^2 = cos n((p + ip)cri + smn(ip + il>)&2, v-3 = 03, 0i 
being the Pauli matrices and n the winding number of the solutions, n — 1, 2, .... In the flat space 
limit, the reduced action density describes a U(l) Higgs like model with two effective Higgs fields 
X A and £ A (A — 1, 2), coupled minimally to the U(l) gauge connection (a r ,ag) [3"T] . 

To remove the U(X) residual gauge freedom of the connection, we impose the usual gauge 
condition d r a r + -dgag = . 
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2.2 Boundary conditions 

In this paper we shall consider both globally regular and black hole solutions of the field equations 
©, (HJ)- The boundary conditions satisfied at infinity and at 9 — 0, ir/2 is the same in both cases, 
and are found from the requirements of finite energy and regularity of solutions. At r — > oo one 
imposes 



a r = Q,a e = -2m, X A = (-l)^n Z'Z , f A = -» ( ) , / = I = p = - = 1 , (8) 



sin 2to# \ c a I sin 2to# 

cos2m9 J \ cos2m6 

with m a positive integer. The following boundary conditions holds for gauge potentials at 9 = 

a r = -d r e, a e = -d e e, x' = 0, ^=0, d eX 2 = 0, £ 2 = -n , (9) 
n n 

while for 6* = ir/2 one imposes 

a r = -£> rX \ ae = -d eX 1 , X 1 = 0, ^ = 0, x 2 = -n, ^ 2 = 0. (10) 
n n 

The boundary conditions for the metric functions at 9 = Q are <9#/ = = = deV = 0, 
and agree with the boundary conditions on the 9 = ir/2 axis. There are also elementary flatness 
requirements which imposes for the metric functions s = I at 9 = and s = p at 9 = n/2. 

To obtain globally regular EYM solutions with finite energy density we impose at the origin 
(r = 0) the boundary conditions 

.i _ / \ tA ( 



a r = , ag = , X A = _ n J , r = ( _ n ] , ft-/ = fts = ft^ - d r p - 0. (11) 

The black hole configurations possess an event horizon located at some constant value of the radial 
coordinate > 0, where the following boundary conditions are imposed 

a r = , d r a e = , drX A = , d r £, A = , f = s = l=p = 0. (12) 

For m = n = 1, these are the spherically symmetric solutions discussed in [3T] , [TO] , [22] . In this 
case the metric functions present no angular dependence, with I = p = s, while ag = w(r) — 1, 
a r = 0, x 1 = -£* = - 1) sin26», \ 2 = -{w{r) - 1) cos 2 9-1, <f = -(u>(r) - 1) sin 2 0-1. 

2.3 Physical quantities 

The mass M of solutions is the conserved charge associated with the Killing vector v — d/dt and 
can be read from the asymptotic expression of the <7 tt -component of the metric tensor 

-fe=/=i-i£+o£). w 

The mass can also be expressed as an integral |32j over the 3-sphere at spacelike infinity, 



bi£^ E - (14) 



M = 

16 

The topological charge of the particle-like solutions as evaluated in [31] is 

<Z^[l-(-l)"> 2 , (15) 

such that the Pontryagin charge is nonzero only for odd to, being equal to n 2 . For even values of 
to, the solutions will describe soliton-antisoliton bound states. 
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To evaluate the Hawking temperature and entropy of the black hole solutions, we use the 
following expansions of the metric functions at the horizon 

The zeroth law of black hole physics states that the surface gravity k is constant at the horizon of 



the black hole solutions, where k 2 — — (1 / A)g u g li ' {digtt){dj9tt) 



Since from general arguments 



the Hawking temperature Tjj is proportional to the surface gravity k, Tjj = ^/(27r), we obtain 
the relation 

= W) 

One can show, with help of the (r (9)-component of the Einstein equations which implies /2S2.0 = 
252/2,6, that the temperature Th, as given in (|16| . is indeed constant. 
For the line element © , the area A of the event horizon is given by 



A = A^ h r /2 d B^eco B eM ^^ e K (17) 
Jo V -' 2 ( ' 

According to the usual thermodynamic arguments, the entropy S is proportional to the area A, 
S = A/AG. 

We mention here also the Smarr-type relation which follows from (|14l) together with Einstein 
equations 

2 47r 2 f°° r /2 x/Tvs 1 

-M = T H S~— dr dQ suiQ cos 0-^— j— (Tl — -T). (18) 



3 6 J rh J / v 1 3 

This relation has been used in practice the verify the accuracy of the numerical computation 



3 Properties of the solutions 

The numerical calculations in this paper were performed by using the software package CADSOL, 
based on the Newton- Raphson method [33J- In this approach, the field equations are first discre- 
tised on a nonequidistant grid and the resulting system is solved iteratively until convergence is 
achieved. In this scheme, a new radial variable x = r/(l + r) (or x = 1 — ru/f for black hole 
solutions) is introduced which maps the semi-infinite region [0, 00) (or [r/j, 00)) to the closed region 
[0,1]. 

For any set of boundary conditions, we have found that the numerical iteration fails to converge 
for T2 = 0. Thus, similar to the spherically symmetric case, no reasonable EYM-p = 1 solutions 
with bi-azimuthal symmetry is likely to exist. This agrees with the physical intuition based 
on a heuristic Derick-type scaling argument (although a rigorous proof exists for the spherically 
symmetric limit only [H1H2]). ^ ^ s * ne -P = 2 YM term, scaling as L~ 8 , which enables the existence 
of configurations with finite mass and well defined asymptotics. 

As in the spherically symmetric case |10j . dimensionless quantities in this model are obtained 
by rescaling the radial coordinate r — > (T2/ri) 1 / 4 r. This reveals the existence of one fundamental 
parameter which gives the strength of the gravitational interaction a 2 = t± 2 (16ttG /t% )• Thus 
without loss of generality, one can fix the values of n and t 2 to some arbitrary positive values and 
construct the solutions in terms of a. We use this property to set in the numerical computation 
r\ = T2 = 1 for m = 1 solutions and t\ = 1, t-i = 1/3 for m = 2 configurations. 

For any set (to, n), the limit a — * can be approached in two ways and two different branches 
of solutions may exist. The first limit corresponds to a pure p = 1 YM theory in a flat background 
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Figure 1: The mass M divided by n 2 is shown as a function of a for globally regular EYM solutions: (a) 
for m = 1, and (b) for m — 2. 



(i.e. no gravity and no p — 2 YM terms), the solutions here replicating the (multi-)instantons 
and composite instanton-antiinstanton bound states discussed in [31J . The other possibility cor- 
responds to a finite value of G as T\ — > 0. Thus, the second limiting configuration is a solution of 
the truncated system consisting of p = 2 YM interacting with gravity, with no p = 1 YM term. 

3.1 Particle- like solutions 
3.1.1 in — 1 configurations 

The m = 1 configurations carry a topological charge n 2 and describe (multi-)solitons. The n = 
1 spherically symmetric case was discussed in [1(5] in a Schwarzschild coordinate system. We 
repeated the numerical analysis of [10] using the isotropic coordinate system (J6)). In the spherically 
symmetric limit only two of the functions in ^ are independent, / and s = I = p. The dominant 
term at the gravity decoupling limit a — > is the F(2) term, the YM solution being the well known 
BPST instanton [33]. When a increases, these solutions get deformed by gravity and the mass M 
decreases (see Figure la). At the same time, the values of the metric functions / and s at the origin 
decrease, as indicated in Figure 2. This branch of solutions exists up to a maximal value a max 
of the parameter a. Another branch of solutions is found on the interval a £ [a cr ^, a max \. On 
this second branch of solutions, both /(0) and s(0) continue to decrease but stay finite. However, 
a third branch of solutions exists for a € [a: cr m, <x C r(2)] > on which the two quantities decrease 
further. A fourth branch of solutions has also been found, with a corresponding a cr (3) close to 
a cr (2). Along this succession of branches, the values of the metric functions / and s at the origin 
continue to decrease. 

On the other hand, the mass parameters do not increase significantly along these secondary 
branches. This behaviour with respect to the parameter a is the same as that which was found 
in [TU], for the metric function a(r) at r = 0. An analytic explanation of these results was given 
in [T2] , where the observed oscillatory behaviour of these functions at r = was characterised as 
a conical fixed point. 

The n > 1 non-spherically symmetric solutions are constructed by starting with the known 
spherically symmetric configuration and increasing the winding number n in small steps. The 
iterations converge, and repeating the procedure one obtains in this way solutions for arbitrary 
n. The physical values of n are integers. We have studied m = 1 solutions with n — 2,3. As 
expected, the general features of the spherically symmetric solutions are the same for all n > 1 
multi-solitons. Like for the Yang-Mills dilaton (YMd) model discussed in [29], when a is increased 
from zero, a branch of gravitating solutions with winding number n emerges smoothly from the 
corresponding F(2) flat space multi-instanton solution. 
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Figure 2: The value at the origin of the metric functions / and s are shown as a function of a for m = 1 
particle-like solutions with n = 1, 2. 



This branch extends up to a maximal value a max (n) of the coupling constant a, beyond 
which the numerical iteration fails to converge. The value of a ma x(n > 1) is smaller than the 
corresponding value in the spherically symmetric case. For example, we find numerically a max {n = 
2) sa 0.412 while the corresponding value for n = 1 is a max w 0.571. For all values n > 1 
we considered, the limiting solutions at a max (n) has no special features. A secondary branch, 
extending backward in a emerges at a max (n). However, the numerical accuracy deteriorates 
drastically for the secondary branch of solutions around some critical value a cr ~ 0.38. Our 
numerical results in this case are less conclusive, the properties of these configurations requiring 
further work. We notice, however, that the value at the origin of all metric functions decreases 
along these branches, as seen in Figure 2 0. We expect that the oscillatory pattern of <?tt(0) 
arising from the conical fixed point observed for the spherically symmetric m = 1, n = 1 solutions, 
will also be discovered for the n > 1 solutions here. However, the construction of the secondary 
branches of solutions is a difficult numerical problem beyond the scope of the present work. 

In all cases we have studied, the metric functions /, I, p, s are completely regular and show no 
sign of an apparent horizon, while I and p have rather similar shapes. The angular dependence of 
the metric functions is rather small, although it increases somewhat with n. The gauge functions 
a r , ag 7 x A > £ A look very similar to those of the YMd solutions presented in [29]. Both \x\ = 
((X 1 ) 2 + (X 2 )) 1/2 and |£| = ((C 1 ) 2 + (£ 2 )) 1/2 possess one node on the 9 = and 6 = vr/2 axis, 
respectively. The positions of these nodes move inward along the branches. 

It is also interesting to note that for the m = 1 solutions, the mass per unit charge of the 
gravitating multisoliton solutions is lower than the mass of a single particle, see Figure la. Thus 
these multisolitons are gravitationally bound states. This case resembles the situation found for 
d = 4 gravitating EYMH monopoles with a vanishing or small Higgs selfcoupling [35] . 

3.1.2 nri == 2 configurations 

The m — 2 configurations reside in the topologically trivial sector. These solutions can be 
thought of as composite systems consisting of two components which are pseudoparticles of Chern- 
Pontryagin topological charges ±n 2 . This type of solutions have no spherically symmetric limit. 
The position of each constituent can be identified according to the location of the maxima of the 
energy density. Also, the structure and location of the nodes of the (effective Higgs) scalar fields 
nicely reveal the evolution and the types of the solutions present at the respective values of the 
gravitational strength. 

2 Note that the values at the origin of all metric functions exhibited in this paper correspond to f(r = 0,8 = 
0), s(r = 0, 8 = 0). This restriction is reasonable since for all solutions with bi-azimuthal symmetry that we have 
found, the metric functions at r = present almost no dependence on the angle 8. 
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Figure 3: The same as Figure 2 for m = 2 particle-like solutions with n = 1, 2. 



As in the case of the m = 1 configurations, coupling with gravity yields various branches of 
gravitating solutions which, however, have different limits depending on the values of the topolog- 
ical charge n 2 of the constituents. Also, their behaviours as functions of the gravitational coupling 
a differ from those with m = 1 presented above. 

n = 1 

There is a certain similarity between the properties of the 4+1 dimensional YMd model studied 
in [29] . and the model under consideration here. As in the former case, we find that in the limit 
a — > resulting from G — > 0, no solution with n = 1 exists, i.e. that in the gravity decoupling 
limit no such solution exists. On the other hand, we know from the work of [TU] that in the flat 
space limit the EYM solution of this model reduces to the BPST instanton [34] of the p = 1 
(usual) YM model, so that in this limit the p = 1 YM term dominates over the p = 2 term. Thus 
the nonexistence of a m = 2 , n = 1 solution here in the gravity decoupling limit implies that 
there should exist no such solution in the 4 + dimensional p = 1 YM model on flat space. This 
is precisely what was found in [31] . 

In the other limit of a — > however, when both r% — > and the gravitational coupling G remain 
finite, such solutions exist. It turns out that in this limit, it is the p = 2 term which dominates 
over the p = 1 YM term. The characteristic feature of this this configuration is that both nodes 
of the effective Higgs fields |x| and |£| merge on the 9 = 7r/4 hypersurface. From this limiting 
configuration, a branch evolves as a increases. Along this branch the nodes move towards the 
symmetry axes, p and a, respectively (with p — rsin#, a — pcos9), forming two identical vortex 
rings whose radii slowly decrease, while the separation of both rings from the origin also decreases. 
The evolution of the solution along this branch can be associated with the increase of the coupling 
n, while T2 and the gravitational coupling G remain fixed. This reproduces the corresponding 
pattern in the YMd system [29]. Note that there is a difference between the evolutions of the 
configurations we are considering here in this 4+1 dimensional theory, and the behaviour of 
the gravitating multimonopoles or the monopole-antimonopole solutions of the gravitating YMH 
system 3 + 1 theory [35l [36] . Although the latter also feature different branches, the evolution 
along those branches is usually associated with the increasing of the gravitational coupling G on 
the lower mass branch, and, the decreasing of the VEV of the Higgs field on the upper mass 
branch. More importantly, the m = 2 , n — 1 solution in that case does have a gravity decoupling 
limit. Thus, the gravitating solutions of the 3 + 1 YMH theory usually are linked to flat space 
configurations, while the solutions discussed here clearly do not have a flat space limit. 

On the p — 2 branch (where the p = 2 term ^Fmnrs dominates) of five dimensional EYM 
m = 2 , n = 1 solutions, the gauge functions a ri ag as well the metric functions / and s are almost 
^-independent, whereas the metric functions / and p possess reflection symmetry with respect to 



8 




0.675 0.7 0.725 0.75 



m=2 n=l 



(2) 



\ m=2 n=2 



(1) 



0.2 0.3 0.4 0.5 

a 



Figure 4: The position of nodes (Zq , Zg 2 ') i s presented for m — 2 particle-like solutions with n = 1, 2. 



= 7r/4 axis. As seen in Figures 1 and 2, the mass of the gravitating solutions on this branch 
decreases, as well as the values at the origin of the metric functions. 

At the critical value a ~ 0.672, the node structure of the configuration changes and both vortex 
rings shrink to zero size, two isolated nodes appearing on each symmetry axis. This transition 
means that the p = 1 term T\ lN becomes dominant. This secondary branch has a small extension 
in a up to the maximal value a max ~ 0.6765, beyond which we could not find regular gravitating 
solutions. We found instead that this branch merges here with the second, p = 1 branch, which 
evolves backwards in a as the value of the metric function /(0) continues to decrease. 

The evolution along this short branch can be associated with the decrease of the coupling 
constant T2 relative to n, as the gravitational coupling G remains fixed. For this branch the 
relative distance between the nodes increases, one lump slowly moving towards the origin and the 
other one moving in the opposite direction. This branch persists up to a value of the coupling 
constant a cr ~ 0.6665, where a critical solution is approached. Due to severe numerical difficulties 
encountered here, we could not clarify the properties of this critical solution further. As a — * a cr , 
the metric function /(0) takes a very small value, f(0) ~ 10~ 3 , while s(0) remains one order of 
magnitude larger (see Figure 3) . At the same time, the Lagrangian density and the mass of the 
configuration remain finite at that point. The critical behaviour observed here resembles the case 
of the gravitating 4+1 EYM vortices in the model consisting only of the p = 1 YM term [5T] . It is 
tempting to speculate that, similar to case in [21] . the solution splits into two parts: a non-singular 
interior region with a special geometry (so-called throat) and an exterior asymptotically flat region 
where two pseudoparticles are located. However, another parametrisation of the metric, differing 
from ([6]) (and possibly even a different numerical approach) appears to be necessary to clarify 
these aspects. 

n = 2 

This configuration also resides in the topologically trivial sector and can be considered as consisting 
of two pseudoparticles of charges ±2 2 . In this case the interaction between the non Abelian 
matter fields becomes stronger than in the case of ±1 constituents, resulting in a different pattern 
of possible branches of solutions. Indeed, as in the case of the 4+1 dimensional YMd system 
[29] . we observe two different branches of gravitating solutions, both linked to the a — * limit. 
The lower branch, on which the p — 1 YM term dominates, emerges from the corresponding flat 
space solution of the pure YM theory with vanishing p = 2 term. Varying a along this branch is 
associated with the decrease of t\ , at fixed t 2 and fixed gravitational coupling G. 

For small values of a the corresponding m = 2 , n = 2 solutions possess two (double) nodes of 
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Figure 5: The energy density e = — T\ is shown as a function of the coordinates p — rsin#, a — pcosO 
for a m — 1, n = 2 EYM black hole solutions with a = 0.2, rn = 0.5. 

the fields \x\ and |£| on the p and a symmetry axes, respectively. The locations of nodes correspond 
to the locations of the two individual constituents and the action density distribution possesses 
two distinct maxima on the 9 = 7r/4 axis. As a increases the mass of the solution increases and 
both pseudoparticles move from spatial infinity towards the origin. For values of a smaller than 
a cr ~ 0.635 along this branch, the energy of interaction between the individual pseudoparticles 
is relatively small and both constituents remain individual. We observe that, as the coupling 
constant approaches this critical value from below, the energy of interaction rapidly increases and 
both pseudoparticles form a bound state, as seen from Figures 3, 4. 

This branch extends further up to a maximal value a ma x — 0.7265 where it bifurcates with 
an upper p — 2 branch which extends all the way back to a = 0. Varying a along this branch is 
associated with the increase of ti relative to t\, as gravitational constant G remains fixed. 

Along the upper branch, as a slightly decreases below a max , the inner node inverts direction of 
its movement toward the outer node which still moves inwards. Thus, both nodes on the symmetry 
axis rapidly approach each other and merge forming a two vortex ring solution at a ~ 0.708. The 
action density then has a single maximum on = 7r/4 axis. As a decreases further both nodes 
move away from the symmetry axis and their positions do not coincide with the location of the 
maximum of the action density. Further decreasing a results in the increase of the radii of the 
two rings around the symmetry axis, and in the limit a — > the rings touch each other on the 
6 = 7r/4 hyperplane. 

3.2 Black hole solutions 

According to the standard arguments, one can expect black hole generalisations of the regular 
configurations to exist at least for small values of the horizon radius r^. This is confirmed by the 
numerical analysis for m = 1, n = 2. Several black hole solutions with m = 2, n = 1 have been 
also constructed, with a lower numerical accuracy, however. 

As discussed in |10j spherically symmetric m = 1, n = 1 black hole counterparts exist for any 
regular solution with the same amount of symmetry. Starting for a given cto < a max from a = 
first branch regular solution, one finds a branch of black hole solutions extending up to a maximal 
value of the event horizon radius = r™ ax . When increases, both the mass and the Hawking 
temperature increase. The value of r h ^ max -^ depends on a. The Hawking temperature decreases on 
this branch, while the mass parameter increases; however, the variation of mass is relatively small. 
The corresponding picture for secondary branches is more complicated and will not be discussed 
here. 
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The numerical construction of nonspherically symmetric black hole solutions appears to be 
more difficult than in the globally regular case. However, our numerical results indicate that the 
m = 1, n > 1 black hole solutions with bi-azimuthal symmetry follow this general pattern. First, 
black hole solutions seem to exist for all values of a or which regular configurations could be 
constructed (here we restrict again to first branch solutions). Also, it appears that black hole 
solutions exist only for a limited region of the (r^,a) space. However, for a given value of a, it 
is very difficult to find an accurate value of r™ ax . An approach to this problem with a different 
method appears to be necessary. 

These solutions possess a regular deformed S 3 horizon. The energy density has a pronounced 
angle-dependence, with a maximum on the 9 = 7r/2 hypersurface. Figure 5 shows a three dimen- 
sional plot of the energy density of a m — 1 , n = 2 black hole with a — 0.2, = 0.5 as a function 
of the coordinates p — rsin#, a — rcosO. With increasing the winding number n, the absolute 
maximum of the energy density residing on the p = a axis, shifts inward. The metric and gauge 
functions possess a nontrivial angular dependence at the horizon. 

Outside their event horizon, these black holes possess nontrivial non Abelian fields. Therefore 
they represent a further counterexample to the d = 5 no-hair conjecture. Also, these bi-azimuthally 
symmetric black holes clearly show that the higher dimensional static black hole solutions need 
not be spherically symmetric. 

4 Conclusions 

Motivated by the recent interest in gravitating solutions in higher dimensional spacetime, we have 
studied static, bi-azimuthally symmetric solutions with non Abelian fields in d = 4 + 1 spacetime 
dimensions. Our solutions are akin to the static, axially symmetric EYM configurations in d = 4, 
studied exhaustively in [37], [351 135] . Our choice of bi-azimuthal symmetry is motivated by our 
desire to reduce the boundary value problem to a two dimensional one. An alternative symmetry 
imposition resulting in a two dimensional residual system would be imposition of SO (3) spherical 
symmetry in the 3 dimensional spacelike dimensions like in 28 . We have eschewed this alternative 
for purely technical reasons (see footnote 1 ). 

The regular and black hole solutions presented are natural generalisations of the known [TU] 
d = 5 EYM spherically symmetric globally regular and black hole solutions. Like the former they 
are asymptotically flat, finite mass solutions, that describe nontrivial gravitating magnetic gauge 
field configurations. Our d = 5 EYM configurations are the first d > 5 dimensional static solutions 
in the literature, which are not spherically symmetric. 

In the case of particle like solutions, which we have studied much more intensively than their 
black hole counterparts, their dependence on the effective gravity coupling a is analysed nu- 
merically in some detail. By and large this is qualitatively very similar to that for the YMd 
solutions [29j in 4 + 1 dimensions, except that here we have four metric functions to keep track of, 
as opposed to the single dilaton field in the previous case [29]. We have studied regular solutions 
with m = 2 , n = 1 and m = 2 , n = 2 in detail, numerically. 

Just as in the YMd case, here too there exists am = 2,n = l solution on the branch where the 
p = 2 YM term dominates, while on the other branch, where the p = 1 YM term dominates, such 
a solution is absent. As it turns out the p = 1 YM term dominates in the gravity decoupling limit, 
which is consistent with our knowledge that this model in 4 + dimensions does not support 31J 
a m = 2 , n = 1 solution. 

Another qualitative feature of 5 dimensional EYM solutions that is confirmed here is the 
occurrence of a conical singular behaviour with respect to the dependence of the metric functions 
on a. This features the oscillatory picture first discovered for the m = 1, n = 1 spherically 
symmetric solutions in |10j and analysed in [12j . which are found also here for the m — 1 , n > 1 
case. 

As compared to the d = 4 case [351 I3H], we expect the existence of a much richer set of 
nonspherically symmetric EYM solutions in d — 5. The configurations studied here represent only 
the simplest, asymptotically flat type of d = 5 nonspherically symmetric gravitating nonabelian 
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solutions. For example, it is known that d — 5 Einstein gravity coupled to Abelian fields presents 
black ring [ID] solutions. These solutions have an horizon topology S 2 x S 1 and approach at 
infinity the flat A4 5 background, as is the case with our solutions. It would be interesting to 
construct non Abelian versions of the U(l) black ring solutions. A black ring can be constructed 
in a heuristic way by taking a black string, bending the extra dimension and spinning it along 
the circle direction just enough so that the gravitational attraction is balanced by the centrifugal 
force. In this framework the (putative) nonabelian black ring would behave locally as a boosted 
black string, e.g. that in [26], with very similar charges and fields. The numerical work involved 
in the construction of non spherically symmetric higher dimensional EYM solutions is, however, 
a considerably challenging task. 
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